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ABSTRACT 

We address three questions regarding solar system planets. What determined their number? Why are 
their orbits nearly circular and coplanar? How long did they take to form? 

Runaway accretion in a disk of small bodies resulted in a tiny fraction of the bodies growing much 
larger than all the others. These big bodies dominated the viscous stirring of all bodies. Dynamical 
friction by small bodies cooled the random velocities of the big ones. Random velocities of small bodies 
were cooled by mutual collisions and/or gas drag. Runaway accretion terminated when the orbital 
$H ' separations of the big bodies became as wide as their feeding zones. This was followed by oligarchic 

Oh, growth during which the big bodies maintained similar masses and uniformly spaced semi-major axes. 

As the oligarchs grew, their number density decreased, but their surface mass density increased. We 
depart from standard treatments of planet formation by assuming that as the big bodies got bigger, the 
small ones got smaller as the result of undergoing a collisional fragmentation cascade. It follows that 
oligarchy was a brief stage in solar system evolution. 

When the oligarchs' surface mass density matched that of the small bodies, dynamical friction was no 
longer able to balance viscous stirring, so their velocity dispersion increased to the extent that their orbits 
crossed. This marked the end of oligarchy. What happened next differed in the inner and outer parts 
of the planetary system. In the inner part, where the ratios of the escape velocities from the surfaces 
of the planets to the escape velocities from their orbits are smaller than unity, big bodies collided and 
coalesced after their random velocities became comparable to their escape velocities. In the outer part, 
where these ratios are larger than unity, the random velocities of some of the big bodies continued to 
rise until they were ejected. In both parts, the number density of the big bodies eventually decreased to 
the extent that gravitational interactions among them no longer produced large scale chaos. After that 
• their orbital eccentricities and inclinations were damped by dynamical friction from the remaining small 

bodies. 

The last and longest stage in planet formation was the clean up of small bodies. Our understanding of 
this stage is fraught with uncertainty. The surviving protoplanets cleared wide gaps around their orbits 
^ \ that inhibited their ability to accrete small bodies. Nevertheless, in the inner planet system, all of the 

material in the small bodies ended up inside planets. Small bodies in the outer planet system probably 
could not have been accreted in the age of the solar system. A second generation of planetesimals may 
have formed in the disk of small bodies, either by collisional coagulation or gravitational instability. In 
the outer solar system, bodies of kilometer size or larger would have had their random velocities excited 
until their orbits crossed those of neighboring protoplanets. Ultimately they would have either escaped 
from the sun or become residents of the Oort cloud. An important distinction is that growth of the inner 
planets continued through clean up, whereas assembly of the outer planets was essentially complete by 
the end of oligarchy. These conclusions imply that the surface density of the protoplanetary disk was 
that of the minimum solar mass nebula in the inner planet region but a few times larger in the outer 
planet region. The timescale through clean up was set by the accretion rate at the geometrical cross 
section in the inner planet region, and by the ejection rate at the gravitationally enhanced cross section 
in the outer planet region. It was a few hundred million years in the former and a few billion years in 
the latter. However, since Uranus and Neptune acquired most of their mass by the end of oligarchy, they 
may have formed before Earth! 

A few implications of the above scenario are worth noting. Impacts among protoplanets of comparable 
size were common in the inner planet system but not in the outer. Ejections from the outer planet 
system included several bodies with masses in excess of Earth after oligarchy, and an adequate number 
of kilometer size bodies to populate the Oort comet cloud during clean up. Except at the very end of 
clean up, collisions prevented Uranus and Neptune from ejecting kilometer size objects. Only Jupiter 
and, to a much lesser extent, Saturn were capable of populating the Oort cloud with comets of kilometer 
size. 
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1. INTRODUCTION 

Modern scenarios for planet formation may be broken 
down into several stages. The growth of the smallest grav- 



itationally active bodies, planetesimals, is mired in con- 
troversy (Lissauer 1993; Youdin & Shu 2002). Orderly 
growth by the merging of planetesimals is followed by run- 
away accretion in which a small fraction of the bodies grow 
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much larger than all the others (Safronov 1972; Wether- 
ill & Stewart 1989). When these big bodies are sparse 
enough so that each dominates viscous stirring in its feed- 
ing zone, runaway growth gives way to oligarchic growth 
during which the big bodies grow in lockstep maintaining 
similar masses and uniformly spaced orbits (Kokubo & Ida 
1998). As oligarchs grow, their orbital spacing increases 
and their number decreases. We investigate how oligarchy 
ends and what happens after it does. The plan of our pa- 
per is as follows. We describe the conditions that pertain 
at the end of oligarchy in §2. We show in §3 that at this 
stage dynamical friction from the small bodies is no longer 
able to balance the mutual stirring of the big bodies. §4 
treats the regularization of the orbits of the big bodies and 
the clean up of small bodies. We summarize our findings 
in §5.1. 

A few definitions arc in order. For simplicity, we con- 
sider two classes of bodies, big ones and small ones, each 
composed of material density p. We denote the surface 
mass density, random velocity dispersion, and radius of 
the former by S, v, and R, and of the latter by a, u, and 
s. The distance from and angular velocity about the Sun 
are given by a and Q. In our numerical estimates, we 
set p — 5.5 gem -3 , the density of Earth, at 1 AU and 
p = 1.5 gem -3 , approximately the densities of Uranus 
and Neptune, at 25 AU. For the condensible fraction of 
the protoplanetary nebula, we adopt the surface densities 
a = 7 gem -2 at 1 AU and a = 1.5 gem -2 at 25 AU. The 
former is just that appropriate to the minimum mass solar 
nebula (Hayashi 1981), but the latter is enhanced sixfold 
relative to it. This enhancement is designed to make the 
isolation mass in the outer planet system comparable to 
the masses of Uranus and Neptune. This is necessary, 
since the timescale for the accumulation of the outer plan- 
ets by coagulation of smaller isolation masses would ex- 
ceed the age of the solar system. The particular value of 
six applies if half the mass had accreted into protoplanets 
by the end of oligarchy. Instead, if most of it had, three 
would be the appropriate enhancement factor. In evalu- 
ating expressions containing the planet mass, M p , we use 
IMo for an inner planet and 15M ffl for an outer one, where 
M® w 6.0 x 10 27 g is Earth's mass. 

It proves convenient to employ a symbol a for the ratio 
between the radius of a body and that of its Hill sphere, 
Rh- 
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234" 1 for a = 1 AU 
K 4ir pa 3 j L 3800 -1 for a = 25 AU 

Note that a is approximately the angle subtended by the 
sun at a. We also make use of the Hill velocity of the 
big bodies, v H ~ £IRh] v h ~ a 1 / 2 w osc , where 

^esc IS the 

escape speed from the surface of a big body. 



proportional to s. We assume that as the big bodies grow 
and viscous stirring intensifies, the small bodies are colli- 
sionally fragmented. Fragmentation lowers their random 
velocities making them easier to accrete. For simplicity, in 
§2.1 we scale all quantities by the values they would have 
if u were fixed at the boundary between the shear and 
dispersion dominated regimes, that is u/vh ~ 1- General- 
ization to other values of u/vh is given in §2.2 along with 
the dependence oiu/vn on s. 

2.1. Conditions at Isolation for u « vh 

With u = vh throughout oligarchy, v < u with equality 
obtaining at isolation (see eqns. [11] & [13] in §3). Thus 
each oligarch accretes material from within an annulus of 
half- width 2.5Rh (Petit & Henon 1986), so the isolation 
radius 
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1.3xl0 3 km fora=lAU 
2.3xl0 4 km fora = 25AU 



Equivalently, the isolation mass 
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The number of oligarchs per unit logarithmic semimajor 
axis is 
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for a = 25 AU 
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Equations (2)-(4) apply for u < Vh- However, from here 
to the end of §2.1 we specialize to u = vh- 

With u = vh, the ratio of the accretion cross section to 
the geometric one is (v esc /u) 2 w a . Thus the timescale 
from the start of oligarchic growth until isolation is 



pa\ V2 f i 5 yr for a = 1 AU 
cr J ~ \ 10 7 yr for a = 25 AU 



(5) 

At isolation, the escape velocity from the surface of an 
oligarch is given by 
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These values are to be compared with the escape veloc- 
ity from solar orbit, 44kms -1 at 1 AU and 8.8kms _1 at 
25 AU. 

Viscous stirring of small bodies by oligarchs at isolation 
results in collisions at speeds w co i <~ Vh, so 
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2. AT ISOLATION 

Here we assume that oligarchy ends when approximately 
half the original surface density has been accreted. 1 We 
refer to this as the epoch of isolation. 

The size, s, and velocity dispersion, u, of the small bod- 
ies are uncertain. They are also closely related: u is set by 
an equilibrium between viscous stirring by the big bodies 
and damping in mutual collisions whose rate is inversely 

1 This assumption is validated in §3 



2.2. Sizes and Velocities of Small Bodies 

The collision rate is directly proportional to the size s 
of the small bodies. To maintain u/vh ~ 1 at isolation 
requires the effective radius of the small bodies to take on 
the particular value 
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10 m for a = 1 AU 
m for a = 25 AU 
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In the shear dominated regime, s < Sb, u/vh ~ s/sb- 
This does not affect the 2.5 Rh half width of an oligarch's 
feeding zone. So i?i so , M- lso , and iVi so are still given by 
equations (2), (3), and (4). Moreover, it co i remains com- 
parable to vh (eq. [7]), the typical random velocity at 
which a small body exits an oligarch's Hill sphere. How- 
ever, as a consequence of the reduced thickness of the disk 
of small bodies, U so oc s provided a 1 / 2 Sfc < s < Sb- There 
is no further reduction of t{ so for s < a 1 / 2 Sfc. We note that 
t- lso can be remarkably small. For s = a 1 / 2 s&, 
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3.1. Shear Dominated Accretion, u < vr 

In this regime the big bodies heat each other and are 
cooled by dynamical friction from the small bodies at the 
rates 2 

1 dv SO _ 9 vh aQ _ 9 

-— ~ a a 10 

v dt pR v pR y ' 

At equilibrium, dv/dt = 0, 

v ~ — vh ■ (11) 
ct 

This justifies the use of the heating rate appropriate to 
v < vh in equation (10). 

3.2. Dispersion Dominated Accretion, u > Vh 

It suffices to consider the regime in which v > Vh as 
well as ii > vh- Under these conditions the big bodies 
heat each other and are cooled by dynamical friction from 
the small bodies at the rates 

(12) 



1 dv Sfi _ 2 fvn\ 4 o-Q _ 2 (VhY 
v dt pR V v J pR V u J 



Equilibrium occurs at 
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3.3. Instability Of Protoplanet's Velocity Dispersion 

As soon as S > ct, the velocity dispersion of the big 
bodies destabilizes. This occurs because the typical rela- 
tive velocity between a big and small body is v > u, so 
equations (10) or (12) are modified to 3 



v dt 



(S — ct) Q, _ 2 (vh > 

-^R— a It 



which follows from multiplying the values in equation (5) 
by a 1 ' 2 . 

In the dispersion dominated regime, s > Sb, u/vh ~ 
(s/sh) 2 / 9 . Thus u co i ~ u, and the width of an oligarch's 
feeding zone u/fl ~ (s/sb) 2 ^ 9 RH- Consequently, to obtain 
values for i?i so , M; so , iVi so , and ti so appropriate to the dis- 
persion dominated regime, we must multiply those given 
in equations (2), (3), (4), and (5) by factors of (s/s^) 1 / 9 , 
(s/sfc) 1,/3 , (s/s b )~ 1/3 , and (s/s h ) 5 / 9 , respectively. 

2.3. Summary 

There are several messages to take away from this sec- 
tion. 

• A short isolation timescale requires accretion of 
small bodies. 

• As the result of viscous stirring by oligarchs, small 
bodies suffer collisions at velocities > Vh that we 
assume are sufficient to fragment them. 

• We note that ti so as used by us measures the du- 
ration of oligarchy. It is quite possible that it is 
shorter than the timescales for orderly and runaway 
growth, the two stages that precede it. 

• Models for Uranus and Neptune imply that each 
planet contains a few M® of hydrogen and helium 
(Guillot 1999). This is consistent with a formation 
timescale ~ 10 7 yr which requires that these plan- 
ets grew by accreting mainly meter size or smaller 
bodies. 

• If collisional fragmentation continues to small 
enough sizes, the disk of small bodies would be op- 
tically thick. Then it would be described by fluid 
rather than by particle dynamics. 

• The isolation mass for the minimum mass solar neb- 
ula is smaller than the planet mass, and by a much 
greater margin in the inner planet system than in 
the outer. 

3. BEYOND OLIGARCHY 

We show below that oligarchy ends when S w a. This 
result applies to accretion in both shear and dispersion 
dominated regimes. 

2 The heating rate given in the literature is proportional to {vjj/v) 2 instead of to Vh/v. Goldrcich ct al. (2004) show that the latter is correct. 

3 Provided v < v csc . 



Thus when S > ct, big bodies are heated faster than they 
are cooled. This marks the end of oligarchy. As v in- 
creases, heating and cooling both slow down, but heating 
always dominates cooling. Eventually the orbits of neigh- 
boring big bodies cross. 

Because it is based on approximate rates for viscous 
stirring and dynamical friction, the criterion, £ ~ ct, for 
the onset of velocity instability is also approximate. Our 
choice of six times the minimum mass solar nebula surface 
density in the outer planet region is based on assuming 
that E = ct at isolation. If instead, at the onset of velocity 
instability, the oligarchs contained most of the mass, the 
appropriate enhancement factor would be slightly above 
three. N-body simulations of oligarch dynamics with the 
addition of accurate analytic expressions for dynamical 
friction can resolve this issue. 

The consequence of the instability in the velocity disper- 
sion differs according to which is larger, the escape velocity 
from the surfaces of the planets that ultimately form or the 
escape velocity from their orbits. The ratio of these two 
escape velocities is 

' 0.3 for a = 1 AU & M p = M e . . 

2.3 for a = 25 AU & M v = 15M ffi " [lb > 
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Before we proceed to discuss these two cases, we stress 
an essential point that is central to the outcome of each. N- 
body planet systems can possess long term stability. This 
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behavior lies outside the realm that naive calculations of 
planetary interactions can describe. We propose that in 
both cases, 1Z < 1 and 1Z > 1, the system of big bodies 
evolves such that the surviving planets have close to the 
smallest spacings allowed by long-term stability. 

3.3.1. Inner Solar System, K< 1: Coalescence 

In regions where K< 1, the big bodies' velocity disper- 
sion increases until it becomes of order the escape velocity 
from their surfaces. At this point they begin to collide and 
coalesce. Coalescence slows as the number of big bodies 
decreases and their individual masses increase. 

The timescale for the formation of planet size bodies 
with radius R p whose orbits are separated by of order a is 
just 



t, 



pR p 



' " " [ ~ctQ ' ~ 10 at ° = 1 AU f ° r Rp = R(1 



(16) 

At a separation of order a, mutual interactions no longer 
produce chaotic perturbations. Indeed, detailed N-body 
simulations of terrestrial planet formation by Chambers 
(2001) produce stable systems on a timescale similar to 

^coag • 

What happens to the small bodies while the big ones 
are colliding and coalescing? A significant fraction col- 
lide with and are accreted by big bodies. Additional small 
bodies are created in grazing collisions between big ones 
(Lcinhardt & Richardson 2002). This ensures that a sig- 
nificant residual population of small bodies persists until 
the end of coalescence. 

3.3.2. Outer Solar System, 1Z 3> 1: Ejection 

In regions where 1Z 3> 1, v reaches the orbital speed £la. 
Some fraction of the big bodies become detached from the 
planetary system and either take up residence in the Oort 
cloud or escape from the sun. This continues until mutual 
interactions among the surviving big bodies are no longer 
capable of driving large scale chaos. 

We estimate the ejection timescale as 

°" J ' M ° * - 10 9 yr at a = 25 AU . (17) 
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Shoemaker & Wolfe (1984) and Dones et al. (2004) report 
similar timcscales for the ejection of test particles placed 
on orbits between Uranus and Neptune, the former from a 
crude impulsive treatment of scattering and the latter from 
N-body integrations. A shorter timescale might apply if 
bodies were transferred to and then ejected by Jupiter and 
Saturn. A quantitative estimate of the transfer rate may 
be obtained from equation (28). 

As the random velocity of a big body increases, the rate 
at which it accretes small bodies declines. Thus a substan- 
tial surface density of small bodies is likely to remain after 
most of the big bodies have been ejected. In the follow- 
ing section we argue that most of the mass in these small 
bodies eventually is either injected into the Oort cloud or 
escapes from the sun. 

4. COMPLETION 



4.1. Gap Clearing 

Gaps were not important prior to isolation because the 
radial spacing of big bodies was only a few times larger 
than the widths of their feeding zones. But after the pro- 
toplanets achieved large scale orbital stability, their radial 
spacing was much larger than their Hill radii and wide 
gaps would have formed around their orbits. 

Gap formation is driven by the torque per unit mass 
(Goldreich & Tremaine 1980) 



T p ~ sgn(x) 
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(18) 



that a protoplanet exerts on material at distance x = a—a p 
from its orbit. 4 The gap width increases with time accord- 
ing to 
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Here we consider processes that took place at sufficiently 
late times, > 10 8 yr in the inner planet region and > 10 9 yr 
in the outer, that it seems safe to ignore effects of gas drag. 

4 Equation (18) is obtained by a radial smoothing of the torque which has peaks at mean motion resonances 



where icyr = t/10 9 yr, and we neglect the presence of 
other planets. 

Because the disk's viscosity arises from random motions 
excited by the protoplanet, the width of the gap is indepen- 
dent of the collision rate. Gap edges are sharp or diffuse 
depending upon whether collisions damp the amplitudes 
of epicyclic oscillations excited at conjunctions before or 
after their phases decohere (Borderies et al. 1989). The 
former would allow accretion, albeit inhibited, while the 
latter would shut it off altogether. The appendix provides 
additional details about the excitation of random motions 
and the profiles of gap edges. 



4.2. Orbit Regularization 

Either coagulation or ejection is likely to end with the 
surviving big bodies moving on orbits with eccentricities 
and inclinations of order 1Z ~ 0.3 in the inner planet sys- 
tem and of order unity in the outer planet system. The 
former is seen in N-body simulations of the formation of 
terrestrial planets from a few hundred big bodies, with no 
small bodies present (Chambers 2001). Such orbits do not 
resemble those of solar system planets. In reality, dynami- 
cal friction by the residual small bodies tends to circularize 
and flatten the orbits of the surviving protoplanets. We 
can compare the rate at which dynamical friction reduces v 
to that at which big bodies grow by accreting small ones. 
For v > v csc both rates are based on physical collisions 
and are of the same order. However, for u < v < v csc , 
the rate at which v damps exceeds that at which R grows 
by the factor (v csc /v) 2 for vh < v < v csc , a~ 1 (v/vH) for 
ol x I 2 vh < v < vh, and a -1 / 2 for v < a 1 / 2 ^. These com- 
parisons apply to a planet that either cannot or has yet to 
open a gap around its orbit. 

Dynamical friction continues to act after gap opening. 
Angular momentum and energy are transferred between 
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the planet and the disk of small particles by torques that 
the planet exerts at Lindblad and corotation resonances. 
Ward & Hahn (1998, 2003) used the standard torque for- 
mula (Goldrcich & Trcmainc 1980) and concluded that 
the most potent contributions to the damping of eccen- 
tricity and inclination are due to torques at apsidal and 
nodal resonances. They assessed their contributions to be 
larger, by factors of ~ f2/|ri7 and ~ S7/|f2 np |, than those 
from torques at standard first order corotation and Lind- 
blad resonances. 5 However, this result comes with a num- 
ber of caveats, especially in applications to disks in which 
self-gravity dominates pressure in the dispersion relation 
for apsidal and nodal waves. Ward & Hahn (1998, 2003) 
assume that these waves are excited at apsidal and nodal 
resonances, then propagate away and ultimately damp. 
They further assume that the resonances lie farther from 
the planet than the first wavelengths of the waves. But the 
main excitation of these waves may occur off resonance at 
gap edges, and their long wavelengths suggest that they 
may have more of a standing than a propagating wave 
character (Goldreich & Sari 2003). Each of these features, 
and especially the latter, is likely to reduce the rates of ec- 
centricity and inclination damping, but by amounts that 
arc difficult to reliably estimate. 

4.3. Clean Up 

What was the fate of the residual small bodies that re- 
mained after the protoplanets had settled onto stable or- 
bits? At the end of oligarchy, small bodies and protoplan- 
ets contributed comparably to the overall surface density. 
But today the mass in small bodies is much less than that 
in planets. The asteroid belt contains most of the mass 
not in planets inside the orbit of Jupiter, but it totals 
< 1O _3 M0. Our knowledge of small bodies in the outer 
planet region is less complete, but observations of perihe- 
lion passages of Halley's comet limit the mass of a disk 
at a » 30 AU to be < 10(a/100 AU) 3 M e (Hamid et al. 
1968; Yeomans 1986; Hogg et al. 1991). 

Clean up was both the last and longest stage in solar 
system evolution. It is ongoing in both the asteroid and 
Kuiper belts. The Oort comet cloud was probably popu- 
lated during this stage. We outline our thoughts on clean 
up below. They are speculations based on interweaving 
theory and observation. 

4.3.1. direct accretion of small bodies 

Accretion of small bodies by protoplanets is the most 
obvious mechanism for clean up. The rate at which a pro- 
toplanet gains mass by accreting small bodies with u ~ vh 
from gap edges at \x\ < 2.5i?# is 

4 
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pR p \ Aa J ■ 



(20) 



where Aa is the distance between neighboring planets, <tq 
is the surface density of the small bodies far from the gap, 
and we assume that the gap's surface density profile obeys 
a oc x 4 (see eq.[A-18]). A more relevant expression is that 
for tdcan = aoldao/dt] -1 « 2ir<JoaAa(dM p /dt)~ 1 ; 

2 / a \ 5 

-l 
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clean 



2 x lO" 2 ^ 1 
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a 



Aa 



8 x 10 10 yr 
7 x 10 11 yr 



for a 
for a = 



-- 1 AU 
25 AU 



(21) 



In both the inner and outer solar system, the spacing be- 
tween planets is Aa <~ a/3, so i c i oan ~ 300 Myr for 1 AU, 
and tcioan ~ 3 Gyr for 25 AU. The latter time is uncom- 
fortably long. It would be a factor a 1 / 2 ~ 1/60 smaller 
for u < a x / 2 VH- However, this introduces a new prob- 
lem. Maintaining such a low velocity dispersion requires 
frequent collisions and therefore substantial optical depth. 
This may lead to sharp gap edges and consequently the ab- 
sence of accretion. See the Appendix for more discussion 
of gap structure. 

4.3.2. second generation planetesimal formation 

Towards the end of oligarchy, small bodies attain ran- 
dom speeds of order 10 2 ms _1 (cq. [7]). Collisions at such 
high speeds fragment them to sizes much smaller than a 
kilometer. After orbit rcgularization the protoplanets are 
spaced by many times their Hill radii and viscous stirring 
of the intervening small bodies is considerably weaker. An 
estimate for the rms random velocity of the small bodies 
is given in equation (A-4). With our standard parameters 
it yields 



1/2 



3/2 



0.1ms- 1 
0.3ms- 1 



for a = 1 AU 
for a = 25 AU . 

(22) 

where |x| is radial distance from the protoplanet. Even 
these small rms velocities are likely to be much larger than 
the mean random velocities. That is because the proto- 
planet 's torque is concentrated at discrete mean motion 
resonances, and the nonlinear disturbances it raises damp 
locally (Goldreich & Tremaine 1978). These strongly 
stirred regions near resonances make the dominant con- 
tributions to u rms . A semi-quantitative discussion of this 
point is provided in the appendix. 

Do larger bodies, referred to here as planetesimals, form 
under the conditions that prevail after orbit rcgulariza- 
tion? We are unable to answer this question with confi- 
dence. Instead, we critique the difficulties faced by coagu- 
lation and gravitational instability, the leading candidates 
for planetesimal formation. 

by coagulation 

Without gravitational focusing, coagulation is a lengthy 
process. To double its mass, a body would have to pass 
through the disk a minimum of ps/a times. A potential 
problem is that a small body's rms random velocity esti- 
mated from equation (22) is greatly in excess of the escape 
velocity from its surface, 



^esc ~ a 



0.7(s/km)ms- 1 
0.4(s/ km) ms -1 



for a = 1 AU 
for a = 25 AU . 

(23) 

This would imply that collisions lead to disruption rather 
than to coalescence. Only bodies larger than 



^crit 



Ml 



M Q M d \x\ 3 



5 x 10 2 km 
2 x 10 5 km 



for a = 1 AU 
for a = 25 AU 



(24) 



5 The symbols ro and Q np denote apsidal and nodal precession rates. 
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have w osc > u rms , where in the numerical evaluation, we 
have set the disk mass, Md ~ aa 2 , equal to the planet 
mass, M p , and \x\ = a. 

It might be argued that equation (22) does not apply, 
that chaotic stirring would not occur far from a planet. 
This is certainly true for perturbations from a single planet 
moving on a nearly circular orbit. However, N-body calcu- 
lations by several groups show that stable orbits between 
planets are rare; even those initialized with low eccen- 
tricities and inclinations invariably become planet orbit 
crossers (Gladman & Duncan 1990; Hohnan & Wisdom 
1993; Grazier et al. 1999). Nevertheless, there are a cou- 
ple of reasons to wonder whether coagulation might still 
occur. None of the N-body calculations investigated the 
stability of orbits with initial random velocities as small 
as a few meters per second, and none of them included the 
small amount of damping that passage though the particle 
disk would cause. 

by gravitational instability 

Gravitational instability is another possible mechanism 
for the formation of second generation planctcsimals. It 
has the virtue of being very fast. However, it also faces a 
problem. The formation of solid bodies by gravitational 
instability requires the particle disk to be optically thick. 
Observations of thermal infrared radiation from solar type 
stars constrain the frequency of protoplanetary systems 
with optically thick disks. 

Suppose that the random velocity of the small bodies 
falls below the limit for gravitational instability. That is, 

u < Ustab ~ ■nGa/n 

^ flOcms -1 for a = 1 AU , (25) 
I lms- 1 for a = 25AU . ' 
Gravitational instabilities convert potential energy into ki- 
netic energy of random motions. The development of non- 
linear overdensities requires this energy to be dissipated at 
the collapse rate ~ Otherwise the random velocity dis- 
persion would be maintained near the margin of stability, 
that is u ~ u s tab (Gammic 2001). Inelastic collisions are 
the only option for dissipating energy in a particle disk. 
For the collision rate to match the collapse rate, the parti- 
cle disk would have to be optically thick, <j/(ps) > 1. An 
optically thick particle disk might result from a collisional 
fragmentation cascade. 

The maximum size of a solid body that can form by col- 
lapse without angular momentum loss in a gravitationally 
unstable disk is 

« _ ^-3/2^ f 50 m for a = 1 AU. 

S * p I 2 km for a = 25 AU. 1 ' 

Rapid damping of random velocities suggests that this is 
the size of first bodies that will form by gravitational in- 
stability. Since the escape velocity from their surfaces is 
ttstab; mutual interactions could maintain their random ve- 
locities at an adequate level to stabilize the disk. 

4.3.3. Inner Solar System 

We assume that most of the mass contained in small 
bodies at the end of coalescence ended up in planets and 
that only a small fraction fell into the sun or was ejected 
by Jupiter. This assumption warrants scrutiny, but that 
will not be done in this paper. 

6 For the sake of argument, we assume that the accretion rate would 



The timescale for clean up by the accretion of small bod- 
ies, as given in equation (21), could be comparable to or, 
for u < vh, even shorter than that for coagulation and 
orbit regularization. However, this should not be taken to 
imply that second order planetesimals did not form during 
clean up. 

4.3.4. Outer Solar System 

difficulties with accretion 

Accretion of the small bodies would be the simplest 
solution to clean up. Estimates based on equation (21), 
which assumes u ~ vjj, suggest that it would take a time 
comparable to the age of the solar system for Uranus and 
Neptune to clean up the region between them which has 
Aa ~ a/3, and far longer for Neptune to clean up material 
from outside its orbit where the gap size would be larger 
(see [19]). Although for u < a x l 2 vn the accretion rate 
would be a factor a" 1 / 2 larger, it would require the disk of 
small bodies to maintain a substantial optical depth. This 
might result in sharp gap edges and a negligible accretion 
rate. Given those uncertainties and our crude estimates, 
we cannot exclude the possibility of accretion. 

A more serious issue for our scenario concerns the 
amount of material that might have been accreted after 
isolation. Could Uranus and Neptune have acquired most 
of their mass during clean up? Suppose the initial surface 
density was only twice that of the minimum mass solar 
nebula and that half remained in the form of small bod- 
ies at isolation. Then the isolation mass would have been 
about one tenth the mass of the outer planets. After a 
fraction of the big bodies were ejected, dynamical friction 
from the small ones would have damped the random ve- 
locities of the survivors. These would then have resumed 
accreting small bodies and, once their masses had grown 
sufficiently, their velocity dispersion would have again be- 
come unstable. This cycle would have repeated until all 
the small bodies were accreted. 6 The end result would 
not have been very different from that of our preferred 
scenario. However, unless the original surface density ex- 
ceeded twice that of the minimum mass solar nebula, the 
repeated ejections would have left too little mass to form 
planets as large as Uranus and Neptune. Also, without 
the formation of a second generation of planetesimals the 
connection to comets would be lost (see below). 

conditions for ejection 
Ejection is the alternative to accretion. Our story im- 
plies that up to ~ 100M® of small bodies was ejected in 
connection with the formation of Uranus and Neptune. 
Such a large mass ejection aided by Jupiter and Saturn 
would have been accompanied by a substantial outward 
migration of Uranus and Neptune(Fernandcz & Ip 1984). 
It might even have moved them outside the orbits of most 
of the material from which they formed (Levison & Mor- 
bidelli 2003). 

To examine the conditions needed for ejection, we con- 
sider the fate of a small body with radius s embedded in 
a sea of bodies with radii < s, and with total surface den- 
sity a. It collides with a total mass of order its own on a 
timescale 

^ ~ £ ~ 2 (iTm) fOT 25 AU - (27) 

: been fast enough for this to happen. 
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By comparison, the timescale for a collisionless test par- 
ticle placed on a low eccentricity orbit midway between 
Uranus's and Neptune's orbit to become an orbit crosser 
is 

'"■Vf^"\' 5M yr. (28) 



tc 



2a 



The above Opik-type estimate (Opik 1976) agrees quite 
well with results from N-body simulations (Gladman & 
Duncan 1990; Holman & Wisdom 1993; Grazier et al. 
1999). Hence only bodies with s larger than 



Me 

Mr, 



2 / A \ 5 

/ Aa 



I 2a 



- 2 km 



(29) 



could have become orbit crossers. But ejection takes much 
longer than orbit crossing; t c joct ~ 1 Gyr (eq. [17]). Only 
bodies larger than 



Reject 



0.1^ 

p 



Me 

M n 



500 km 



(30) 



could have been ejected by Uranus and Neptune in the 
presence of a disk of smaller bodies. However, for Jupiter, 
equations (17) and (30) yield i e j oc t ~ 10 5 yr and s c joct ~ 
6 km. 

The Oort cloud is a repository for kilometer size bodies 
that probably formed in and were ejected from the outer 
planet region. Current estimates of the cloud's mass lie 
in the range 1-10 M s (Weissman 1996), with the size of 
pristine comets, 1 km < s < 10 km being a major part 
of the uncertainty. Detailed numerical calculations that 
follow the ejection of test particles from the outer planet 
region show that a few percent end up in the Oort cloud 
(Dones et al. 2004). These, together with the observed flux 
of new comets, are taken to imply that the outer planets 
ejected a few hundred M e of kilometer size bodies. Some 
fraction may have originated in the vicinity of Uranus and 
Neptune and been transferred via Saturn to Jupiter which 
then ejected them. 

Simplified treatments by Shoemaker & Wolfe (1984) and 
Fernandez (1997) as well as N-body simulations by Dones 
et al. (2004) show that 50-80 percent of test particles ini- 
tially placed between Uranus and Neptune are, in fact, 
ejected by Jupiter. But these investigations did not in- 
clude collisional damping whose importance was first rec- 
ognized by Stern & Weissman (2001), and further inves- 
tigated by Charnoz & Morbidelli (2003). When this is 
accounted for, it yields a stronger result: Jupiter and, to 
a much lesser extent, Saturn were responsible for ejecting 
almost all of the kilometer size bodies into the Oort cloud. 
For this scenario to work, kilometer size bodies must have 
formed out of the much smaller collisional debris that ex- 
isted at the end of oligarchy (see §4.3.2). 

5. DISCUSSION 

5.1. Conclusions 

The scenario sketched in this letter addresses some of 
the basic problems in planet formation. 



• The number and orbital spacing of the planets re- 
sulted from an evolution toward stability against 
large scale chaotic perturbations. 

7 The timescale in the outer solar system could have been much shorter if all ejections were done by Jupiter and Saturn 



• After the cessation of chaotic perturbations, dy- 
namical friction by the residual small bodies 
damped the orbital eccentricities and inclinations 
of the surviving protoplanets. 

• Accretion during oligarchy involved small bodies 
created by a collisional fragmentation cascade. This 
stage probably lasted for less than 10 5 yr in the in- 
ner solar system and less than 10 7 yr in the outer 
solar system. 

• The timescale for establishing the final configura- 
tion of planetary orbits was a few hundred million 
years for the inner planet system and a few billion 
years in the outer planet system. It was set by the 
accretion rate at the geometrical cross section in the 
former, and by the ejection rate at the gravitation- 
ally enhanced cross section in the latter. 7 

• Clean up of small bodies is a complicated and 
poorly explored stage of planet formation. Small 
bodies in the inner solar system were incorporated 
into planets. Those in the outer solar system were 
probably ejected by Jupiter and Saturn, but that 
requires a second generation of planctesimal forma- 
tion. 

5.2. Influence of Gas 

We have neglected the influence of gas. Observations of 
young stars indicate that protostellar disks dissipate in a 
few million years (Haisch et al. 2001; Strom et al. 1993). 
We show below that although the presence of gas would 
alter some of our numerical results, it would not affect our 
picture qualitatively. 

Gas drag can provide significant damping for the ran- 
dom velocities of small bodies in addition to that due to 
inelastic collisions. Relative to collisions, it is most effec- 
tive in damping the random velocities of bodies that are 
smaller than the mean free path of the gas molecules. For 
these, its damping rate obeys the same expression as that 
due to collisions, but with the surface density of the small 
bodies replaced by that of the gas. We can account for 
the effects of gas drag by considering s to be an effective 
size for the small bodies that can be less than their true 
size. This is a minor point for our story since, as we have 
emphasized, the true size of the small bodies is highly un- 
certain. Moreover, our main concern is with the stages of 
planet formation that follow velocity instability and these 
probably continue after the gas is gone. 

Rafikov (2003) explored the fast accretion of protoplan- 
etary cores in the presence of gas. His investigation runs 
parallel to the early phases of ours. However, it terminates 
at the onset of velocity instability, when v <~ vr- 

A potentially more significant effect of gas drag that 
was not considered by Rafikov (2003) is its role in damp- 
ing the random velocities of the oligarchs. Ward (1993) 
shows that the gas damping rate can be obtained from 
the damping rate due to small bodies by substituting the 
surface density of gas for that of the small bodies and the 
sound speed of the gas, c s , for the random velocity of the 
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small bodies, it. 8 By stabilizing the oligarchs' random ve- 
locities, gas drag could have enabled them to consume all 
of the small bodies. 

In the inner solar system, it is possible — though highly 
uncertain — that much of the gas survived until isolation. 
Then the full velocity instability of the oligarchs would 
have been delayed until the surface density of gas declined 
to match that contributed by oligarchs. After that the 
oligarchs would have excited their random velocities up to 
their escape speeds. Although most of the small bodies 
would have been accreted before this happened, plenty of 
new ones created in glancing collisions could have damped 
the orbital eccentricities and inclinations of the planets 
that finally formed. 9 

Outer solar system planets, Uranus and Neptune, are 
believed to have collected only a few earth masses of nebu- 
lar gas. So it is likely that most of the gas had disappeared 
prior to isolation in the outer solar system. 

Gas drag must have been more significant in the for- 
mation of Jupiter and Saturn. One might worry that the 
orbital decay of small particles, which are an integral part 
of our scenario, would have been too fast for them to have 
been accreted. Particles with stopping time comparable 
to their orbital time drift fastest. Their orbits decay on a 
timescale f2 _1 (af2/c s ) 2 <~ 10 3 yr. By damping the random 
velocities of small bodies, gas drag can protect them from 
undergoing destructive collisions. This may result in larger 
bodies, for which the drift timcscales are longer. Moreover, 
gas drag could have made the isolation timescale in the 
Jupiter-Saturn region as short as ~ 10 4 yr (see equation 
[9])- 

This research was supported in part by NSF grants 
AST-0098301 and PHY99-07949 and NASA grant NAG5- 
12037. We thank Mike Brown, Eugene Chiang, Luke 
Dones, Martin Duncan, Shigcru Ida, Roman Rafikov, and 
David Stevenson for helpful advice. 

APPENDIX 

In this appendix, we consider how a planet affects a 
disk of small bodies in which it is embedded. We assume 
throughout that the collision time between small bodies 
is longer than the time it takes their epicyclic phases to 
decohere. This applies provided the optical depth of the 
small bodies at a distance \x\ > Rh from the planet is less 
than T cr it ~ (Rh/\x\) 3 . If this condition does not hold, 
i.e., if t > T cr it, then the effective viscosity would be neg- 
ative, and the planet would open a gap with sharp edges 
(Borderies et al. 1989). 

Stirring 

To determine the mean random kinetic energy of the 
small bodies, we balance their energy loss rate in inelastic 
collisions, (1 — e 2 )u 2 /t co \, against the sum of the rates of 
energy gain from direct forcing by the planet, (Q p — fl)T p 
(Borderies et al. 1982), plus viscous dissipation acting on 
the Keplerian shear, v(ad£l/ da) 2 . The symbol e < 1 de- 



notes the coefficient of restitution, assumed to be a de- 
creasing function of impact velocity, fi p is the orbital fre- 
quency of the planet, and T p is the torque per unit mass 
exerted by the planet. For \x\ <C a we obtain 



(1 



^col 



3x „ 
2^> 
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M 2 



where the kinematic viscosity 



(Goldreich & Trcmaine 1978), 10 and 



u n- 
v ~ \l 

^coll 
10 



col 



ps 



(A-l) 



(A-2) 



(A-3) 



In the limit T, p = 0, energy released from the Keplerian 
shear would maintain u at some small value set by the ve- 
locity dependence of e. The velocity dispersion in unper- 
turbed parts of Saturn's rings, u < 1 cms -1 , is a practical 
example. Our interest is in circumstances under which 
forcing by a planet results in an equilibrium value of u 
that is much larger than that produced solely by viscous 
dissipation acting on the Keplerian shear. Under such con- 
ditions 

-^(f) 1/2 feVV (A-4) 



Mr. 



To this point we have proceeded as though all small bod- 
ies were of the same size. However, equation (A-4) applies 
more generally and yields the size dependence of the rms 
random velocity of a small body subject to two limits. It 
must be larger than those which contain most of the mass 
but small enough so that gravitational focusing does not 
enhance its interactions with them. In this generalized 
interpretation, a must be interpreted as the total surface 
density in small bodies of all sizes. 

Resonances 

The planet's torque is concentrated at resonances. 
Equation (A-4) is a spatial average of the torques at mean 
motion resonances. These dominate the heating of the 
random velocities of small bodies. As the rms random ve- 
locity given by equation (A-4) is a spatial average, it is of 
limited utility. 

A more complete picture is obtained by investigating 
the disturbances raised by torques at individual princi- 
pal mean motion resonances (Goldreich & Tremainc 1980). 
Our starting point is the WKB dispersion relation for non- 
axisymmetric waves of angular degree m in a cold, self- 
gravitating disk; 

(w - mQ) 2 = k 2 - 2-xGo\k\ . (A-5) 

For principal mean motion resonances, lu = mtt p , where 
the pattern speed, is equal to the planet's mean or- 
bital angular velocity. At the Lindblad resonance, |A;|=0 
and \x\/a w 2/(3m). 

Density waves are excited at each mean motion reso- 
nance. Their properties have been extensively studied; we 
merely quote a few relevant results (Goldreich & Tremaine 

8 We note that at isolation vh < c s < v CB . Moreover we are assuming that v < c s . 

9 Since it took the inner planets more than 100 Myr to form (eqn. [16]), gas is unlikely to have contributed to regularizing their orbits. 

10 We consider only the case Clt co \ > 1; more generally, Goldreich & Tremaine (1978) show u ~ u 2 t co \ [l + (Slt co i) 2 ] 1 for circular, Keplerian 
rotation. 
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1978). These are specialized to the case of a near Keple- 
rian disk for which k « Cl. A wave propagates away from 
the resonance and the planet at the group speed 

ttG<t 

v = 



Its first wavelength 



Ai_ 

\x\ 



aa 2 a 



1/2 



(A-6) 



(A-7) 



At each encounter with the protoplanet, a disk particle 
receives a kick sufficient to change its orbital eccentricity 

M„ ( a 



Ae 



M G 



© 



(A-8) 



At a Lindblad resonance, successive increments in e sum 
coherently over a time comparable to that during which a 
disturbance propagating at the group velocity crosses the 
first wavelength. The number of encounters that occur in 
this time is 

\x\ 



N 



so that at resonance, 



iVAe 



M 2 



1/2 



M Q aa 2 



(A-9) 



(A-10) 



The nonlinearity of the wave, Act/ct, is of order the ratio of 
the coherent epicyclic excursions to the wavelength. Near 
resonance this gives 



Ai 



M p a 
aa 2 \x\ 



(A-ll) 



In our scenario, the disk and planet have compara- 
ble masses, and both are much smaller than the solar 
mass. Thus the density waves reach order unit nonlinearity 
within their first wavelengths and their first wavelengths 
are much smaller than the distance between neighboring 
resonances. Therefore, the waves only propagate a small 
fraction of the distance between resonances before damp- 
ing. All of the planet's excitation of random velocities is 
concentrated in these narrow regions where they damp. 



where e* is the value of e at which the equilibrium velocity 
dispersion is obtained for stirring by the Keplerian shear in 



the absence of a protoplanet. Since T p 
equation (A-l) also implies that 

,2 
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3QK 



2a(e 



K/x 4 (eq.[18]), 



(A-14) 



In steady-state, the torque per unit mass that the pro- 
toplanet exerts on a small body, K/x , is balanced by the 
viscous torque. Thus 

K 3f2a d{va) 
- . 



(A-15) 



x" 1 2a dx 
where we neglect gradients of 17 and a, since they are much 
smaller than those of va for |a;| <C a. Combining equations 
(A-12), (A-13), (A-14), and (A-15) gives 

d ( (1 - el) o 
dx \{el-e 2 ) \x 



a 



(A-16) 



For circumstances where stirring by the protoplanet is 
much greater than that due to the Keplerian shear, it is 
likely that e <^ e». In this case equation (A-16) yields the 
gap profile 

a cx \x\\ with q = 3 + e 2 J(l - el) . (A-17) 

It is difficult to obtain a reliable estimate for e». Gol- 
dreich & Tremaine (1978) obtained an approximate solu- 
tion of the collisional Boltzmann equation for a model in 
which the particles were represented by smooth spheres 
separated by many times their diameters. They found 
e* « 0.63, which implies 

(joclxl 3 - 66 . (A-18) 

For the purposes of the present paper, we are content to 
approximate this as a oc \x\ 4 . 



Gaps 

A protoplanet clears a gap in the disk of small bodies in 
which it is embedded. Epicyclic motions are excited when 
small bodies pass conjunction with the protoplanet. Pro- 
vided their phases decohere before collisions damp their 
amplitudes, the gap edges will be diffuse rather than sharp. 
Since the accretion rate in the shear dominated limit is 
proportional to the surface density at \x\ ~ 2.5Rh, it is 
important to determine the gap's surface density profile. 
In order to do so we must estimate the order unity coeffi- 
cient relating the kinematic viscosity to the rate at which 
energy per unit mass is dissipated by inelastic collisions. 
We define b through 



M 2 = 



46 u 2 
~9t~ 



y £coi 

Combining equations (A-l) and (A-12), we find 



b=l-el 



(A-12) 



(A-13) 
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